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ABSTRACT 

We summarize calculations of the cross section for top quark production at 



(N 

' hadron colliders within the context of perturbative quantum chromodynamics, 

including resummation of the effects of initial-state soft gluon radiation to all 
• orders in the strong coupling strength. In our approach we resum the universal 

leading-logarithm contributions, and we restrict the calculation to the region of 
■ phase space that is demonstrably perturbative. We compare our approach with 

Oh! other methods. We present predictions of the physical cross section as a function 

of the top quark mass in proton-antiproton reactions at center-of-mass energies 
of 1.8 and 2.0 TeV, and we discuss estimated uncertainties. 

1. Introduction and Motivation 

X 

In this report we summarize calculations carried out in perturbative quantum 
chromodynamics (QCD) of the inclusive cross section for the production of top quark- 
antiquark (tt) pairs in hadron reactionsSii'i. We begin with a discussion of the mo- 
tivation for the inclusion of the effects of intial state soft gluon radiation to all orders 
in the QCD coupling strength, and we review the general formalism of resummation. 
We outline the method and domain of applicability of perturbative resummation that 
we developed in the past yeari'i, and we contrast this approach with other methods^. 
We present predictions of the physical cross section as a function of the top quark 
mass in proton-antiproton reactions at center-of-mass energies of 1.8 and 2.0 TeV. 

In hadron interactions at collider energies, tt pair production proceeds through 
partonic hard-scattering processes involving initial-state light quarks q and gluons 
g. In lowest-order QCD, 0(a%), the two partonic subprocesses are q + q — > t + i 
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and g + g —>■ t + t. Calculations of the cross section through next-to-leading order, 
0(a^), involve gluonic radiative corrections to these lowest-order subprocesses as well 
as contributions from the q + g initial state!. A complete fixed-order calculation at 
order 0(a™), n > 4 does not exist. 

The physical cross section for each production channel is obtained through the 
factorization theorem, 

The square of the total hadronic center-of-mass energy is S, the square of the partonic 
center-of-mass energy is s, and m denotes the top mass. The variable 77 = — 1 
measures the distance from the partonic threshold. The indices ij G {qq, gg} de- 
note the initial parton channel. The partonic cross section &ij(r), m 2 , /1 2 ) is obtained 
commonly from fixed-order QCD calculations!, or, as described here, from calcula- 
tions that go beyond fixed-order perturbation theory through the inclusion of gluon 
resummationMffl The parton flux is ^(y,^ 2 ) = /J ^/^(x, fi 2 )fj/ h2 (y/x, /i 2 ), 
where fi/ hl (x, /x 2 ) is the density of partons of type i in hadron h±. Henceforth, we use 
the notation a(/i) = a s (^)/iT, where /i is the common renormalization/factorization 
scale of the problem. Unless otherwise specified, a = a(// = m). The total physical 
cross section is obtained after incoherent addition of the contributions from the the qq 
and gg production channels. We ignore the small contribution from the qg channel. 

A comparison of the partonic cross section at next-to-leading order with its lowest- 
order value reveals that the ratio becomes very large in the near-threshold region, i.e., 
the ".fT-factor" at the partonic level K(r]) becomes very large as rj — > 0. An illustration 
of this behavior may be seen in Fig. 7 of Ref. [3]. This large ratio casts doubt on the 
reliability of simple fixed-order perturbation theory for physical processes for which 
the near-threshold region in the subenergy variable contributes significantly to the 
physical cross section. Top quark production at the Fermilab Tevatron is one such 
process, because the top mass is relatively large compared to the energy available. 
Other examples include the production of hadronic jets that carry large values of 
transverse momentum and the production of pairs of supersymmetric particles with 
large mass. To obtain more reliable theoretical estimates of the cross section in 
perturbative QCD, it is important first to identify and isolate the terms that provide 
the large next-to-leading order enhancement and then to resum these effects to all 
orders in the strong coupling strength. 

2. Gluon Radiation and Resummation 

The origin of the large threshold enhancement may be traced to initial-state 
gluonic radiative corrections to the lowest-order channels. Consider the subprocess 
+ J (^2) — > t(pi) + t{P2) + g(k). We define the variable z through the partonic 
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invariants^! 



(k 1 + k 2 ) 2 , h = (k 2 -p 2 ) 2 -rn 2 , u\ = (h-p 2 f 



■m 



(l-z)m 2 = s+tt+ut. (2) 



Alternatively, (1 — z) = In the limit that z — > 1, the radiated gluon carries 

zero momentum. After cancellation of soft singularities and factorization of collinear 
singularities in 0(a 3 ), there are left-over integrable logarithmic contributions to the 
cross section associated with initial-state gluon radiation. These contributions are 
proportional to ln(l — z). 

The partonic cross section may be expressed generally as 



<7y(?7, m 2 



dz 



1 



Hij{z,a) 



(3) 



We work in the MS factorization scheme in which the q, q and g densities and the 
next-to-leading order partonic cross sections are defined unambiguously. The lower 
limit of integration, z min = 1 — 4(l + r/) + 4^1 + r/, is set by kinematics. The derivative 
<Jij(j], m 2 , z) = d(alj\r],m 2 , z))/dz, and is the lowest-order 0(a 2 ) partonic cross 
section expressed in terms of inelastic kinematic variables to account for the emitted 
radiation. 

Keeping only the leading logarithmic contributions through 0(a 3 ), we may write 
the total partonic cross section as 



a 



(o+i) 



(rj, m 2 



f 1 dz{l + 



a2C ij ln 2 ( 



(4) 



where C q q = Cp = 4/3 and C gg = Ca = 3. As is illustrated in Fig. 1(a), the 
leading logarithmic contribution, integrated over the near-threshold region 1 > z > 0, 
provides an excellent approximation to the exact full next-to-leading order physical 
cross section as a function of the top quark mass. Although an exact fixed-order 
0(a 4 ) calculation of ti pair production does not exist, we may invoke universality 
with massive lepton-pair production (//), the Drell-Yan process, to generalize Eq. (U) 
to higher order. In the near-threshold region, the kernel becomes 



n 



(0+1+2) 
ij 1 



z, a) 



2aCij£n 2 (l - z) 



2C%ln\\ -z)- -Cijhin^l - z) 



(5) 



The coefficient b 2 = (11C^ — 2n^)/12, and the number of flavors rif = 5. The 
further enhancement of the physical cross section produced by the O (a 4 ) leading 
logarithmic terms in the near-threshold region is shown in Fig. 1. At m = 175 
GeV, we compute the following ratios of the physical cross sections in the leading 
logarithmic approximation: cr}° faff =1.22, and cr-° +1+2 ' ) /o"^ 0+1 ' ) = 1.14. These 
ratios show that the near-threshold logarithms build up cross section in a worrisome 
fashion. The ratios suggest that perturbation theory is not converging to a stable 
prediction of the cross section. 
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Fig. 1. (a) Physical cross sections in the qq channel in the MS scheme. The solid lines denote the 
finite-order partial sums of the universal leading-logarithmic contributions from the explicit 0(a 3 ) 
and 0(a 4 ) calculations for the tt and Drell-Yan cross sections, respectively. Lower solid: a^; middle 
solid: cr( 0+1 ); upper solid: a^ 0+1+2 \ The dashed curve represents the exact next-to-leading order 
calculation for tt production, in excellent agreement with a^ +1 K The dotted curve is our resummed 
prediction, (b) Optimum number of perturbative terms in the exponent with PVR (solid family 
is for PVR, dashed for perturbative approximation, both families increasing, for parametric values 
n = 10,20,30,40). 



The goal of gluon resummation is to sum the series in a n In n (l — z) to all orders 
in a in order to obtain a more defensible prediction. This procedure has been studied 
extensively for the Drell-Yan process^. In essentially all resummation procedures, 
the large logarithmic contributions are exponentiated into a function of the QCD 
running coupling strength, itself evaluated at a variable momentum scale that is a 
measure of the radiated gluon momentum. For example, in the approach of LSvN0, 
the resummed partonic cross section is written as 

*% IRC (V, = r Mmf dze E ^%( V , m\ z), (6) 
where the exponent 

E^z.m 1 ) oc CijOt({l - z) 2/3 m 2 )£n 2 (l - z). (7) 

Different methods of resummation differ in theoretically and phenomenologically 
important respects. Formally, if not explicitly in some approaches, an integral over 
the radiated gluon momentum z must be done over regions in which z —>■ 0. Therefore, 
one significant distinction among methods has to do with how the inevitable "non- 
perturbative" region is handled in each case. Examination of Eqs. (|6]) and ([/]) shows 
that an infrared singularity is encountered in the soft-gluon limit z — > 1: owing to 
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the logarithmic behavior of a(q 2 ), a(q 2 ) oc £n~ l (q 2 / Aq CD ), a((l — z) 2 ^m 2 ) — > oo as 
z —*■ 1. The infrared singularity is a manifestation of non-perturbative physics. In the 
approach of LSvN, this divergence of the integrand at the upper limit of integration 
necessitates introduction of the undetermined infrared cutoff (IRC) /i D in Eq. (|6]), 
with Aqcd < A*o < m - The cutoff prevents the integration over z from reaching the 
Landau pole of the QCD running coupling constant. The presence of an extra scale 
spoils the renormalization group properties of the overall expression. The unfortunate 
dependence of the resummed cross section on this undetermined cutoff is important 
numerically since it appears in an exponential. It is difficult to evaluate theoretical 
uncertainties in a method that depends on an undetermined infrared cutoff. 



3. Perturbative Resummation 

The method of resummation we employlli is based on a perturbative truncation of 
principal-value resummation. The principal-value method (PVR)! has an important 
technical advantage in that it does not depend on arbitrary infrared cutoffs, as all 
Landau-pole singularities are by-passed by a Cauchy principal-value prescription. 
Because extra undetermined scales are absent, the method also permits an evaluation 
of the perturbative regime of applicability of the method, i.e., the region of the gluon 
radiation phase space where perturbation theory should be valid. 

To illustrate how infrared cutoffs are avoided in the PVR method, it is useful 
to begin with an expression in moment (n) space for the exponent that resums the 
£n(l — z) terms!. Factorization and evolution lead directly to exponentiation in 
moment space: 

E(n,m 2 ) = - J dx—^ — ^ j y# [a (Am 2 )] . (8) 

(l-z) 2 

The function g(a) is calculable perturbatively, but the behavior of a(Xm 2 ) leads to a 
divergence of the integrand when Am 2 — > Aq CD . To tame the divergence, a cutoff can 
be introduced in the integral over x or directly in momentum space, in the fashion of 
LSvN. In the principal-value redefinition of resummation, the singularity is avoided 
by replacement of the integral over the real axis x in Eq. (|8|) by an integral in the 
complex plane along a contour P that is symmetric under reflections across the real 
axis: 

E PV (n,m 2 ) = - [d( C -^-± J fg[a(Xm 2 )]. (9) 



The function E PV (n, m 2 ) is finite since the Landau pole singularity is by-passed. In 
Eq. (0), all large soft-gluon threshold contributions are included through the two- loop 
running of a. 
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Equations (§) and (§) have identical perturbative content, but they have different 
non-perturbative content since the infrared region is treated differently in the two 
cases. The non-perturbative content is not a prediction of perturbative QCD. In our 
study of top quark production, we choose to use the exponent only in the region of 
phase space in which the perturbative content dominates. 

We use the attractive finiteness of Eq. @ to derive a perturbative asymptotic 
representation of E(x,a(m)) that is valid in the moment-space interval 

1 < x = Inn < t = — '— . (10) 

This perturbative asymptotic representation is 

N(t)+1 p+l 

E^x, a) ~ E i3 (x, a, N(t)) = 2C i3 £ « P E ■ ( n ) 

p=l j=0 



Here 



_ 6r i ( _l ) ^ 2 P Cp+1 _ j(p _i) ! / j! , (12) 

and T(l + 2) = Z)^ =0 c fc zfc ' where V is the Euler gamma function. The number of 
perturbative terms N(t) in Eq. ( O ^ 



proximation 



£(x,a)-£(x,a,iV(t)) 



is obtainedll by optimizing the asymptotic ap- 
minimum. Because of the range of validity in 

Eq. (0), terms in the exponent of the form a k ln fc n are of order unity, and terms with 
fewer powers of logarithms, a k ln fc_m n, are negligible. This explains why resummation 
is completed in a finite number of steps. In addition, we discard monomials a k In n 
in the exponent because of the restricted leading-logarithm universality between the 
tt and II processes. 

The exponent we use in our calculations of top quark production is the truncation 

N(t)+1 

E lj (x,a,N) = 2C ij J2 a P s P x p+ \ (13) 
P =i 

with the coefficients s p = s p+ i iP = b%~ 2 P '/ 'p(p + 1). The number of perturbative 
terms N(t) is a function of only the top quark mass m. This expression contains no 
factorially-growing (renormalon) terms. The perturbative region of phase space is far 
removed from the part of phase space in which renormalons could be influentiali. 

In Fig. 1(b) we illustrate the validity of the asymptotic approximation for a value 
of t corresponding to m — 175 GeV. Optimization works perfectly, with N(t) = 6, 
and the plot demonstrates the typical breakdown of the asymptotic approximation if 
N were allowed to increase beyond N(t). This is the exponential rise of the infrared 
(IR) renormalons, the (p — 1)! growth in the second term of Eq. (0). As long as n 
is in the interval of Eq. (|T0D, all the members of the family in n are optimized at the 
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same N(t), showing that the optimum number of perturbative terms is a function of 
t, i.e., of m only. 

It is valuable to stress that we can derive the perturbative expressions, Eqs. fllCf), 
(ID]) , and (p!2[), from the unregulated exponent Eq. (||) without the PVR prescription, 
although with less certitude. We discuss this point in some detail in our long paper!. 

After inversion of the Mellin transform from moment space to the physically rel- 
evant momentum space, the resummed partonic cross sections, including all large 
threshold corrections, can be written in the form of Eq. but with Eq. ([5]) replaced 
by 

J ° 3=0 

The leading large threshold corrections are contained in the exponent Eij( cal- 
culable polynomial in x. The functions {Qj(x, a)} arise from the analytical inversion 
of the Mellin transform from moment space to momentum space. These functions are 
produced by the resummation and are expressed in terms of successive derivatives of 
E: Pk(x,a) = d k E(x,a)/k\d k x. Each Qj contains j more powers of a than of x so 
that Eq. ( |14|) embodies a natural power-counting of threshold logarithms. 

The functional form of E^ for tt production is identical to that for 11 production, 
except for the identification of the two separate channels, denoted by the subscript 
ij. However, only the leading threshold corrections are universal. Final-state gluon 
radiation as well as initial-state/final-state interference effects produce sub-leading 
logarithmic contributions that differ for processes with different final states. Among 
all {Qj} in Eq. (0), only the very leading one is universal. This is the linear term in 
Pi contained in Q , that turns out to be Pi itself. Since we intend to resum only the 
universal leading logarithms, we retain only Pi. Hence, Eq. (|i~4|) can be integrated 
explicitly, and the resummed version of Eq. (0) is 

af pert (77, m 2 ) = I*™* dze E ^ H ^ a) a' tJ ( V , m 2 , z). (15) 

The upper limit of integration in Eq. flTSD is set by the boundary between the 
perturbative and non-perturbative regimes. An intuitive definition of the perturbative 
region, where inverse power terms are unimportant, is provided by the inequality 
(t-z)m < 1. This inequality is identical to the expression in moment space, Eq. fllOD , 
with the identification n = In In momentum space, the same condition is realized 
by the constraint that all {Qj}, j > 1 be small compared to Qo- From the explicit 
expressions!! for the {Qj}, one may show that this constraint corresponds to 

Pi (in (-!—), a) <1. (16) 



1 - z, 

Equation ([16]) is equivalent to the requirement that terms that are subleading accord- 
ing to perturbative power-counting are indeed subleading numerically; Eq. (|16|) is the 
essence of perturbation theory in this context. 
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As remarked above, we accept only the perturbative content of principal-value 
resummation, and our cross section is evaluated accordingly. Specifically, we use 
Eq. ( Jl5| ) with the upper limit of integration, z max , calculated from Eq. ([16]). The up- 
shot is an effective threshold cutoff on the integral over the scaled subenergy variable 
rj, but one that is calculable, not arbitrary. Our perturbative resummation probes the 
threshold down to the point 77 > r/ = (1 — z max )/2. Below this value, perturbation 
theory, resummed or otherwise, is not to be trusted. For the top mass m = 175 GeV, 
we determine that the perturbative regime is restricted to i] > 0.007 for the qq channel 
and rj > 0.05 for the gg channel. These numbers may be converted to more readily 
understood values of the subenergy above which we judge our perturbative approach 
is valid: at m = 175 GeV, these are 1.22 GeV above the threshold in the qq channel 
and 8.64 GeV above the threshold in the gg channel. The difference reflects the larger 
color factor in the gg case. A larger color factor makes the non-perturbative region 
larger. (One could attempt to apply Eq. ( |i~5"D all the way to z max = 1, i.e., to r\ = 0, 
but one would then be using a model for non-perturbative effects, the one suggested 
by PVR, below the region justified by perturbation theory.) 

It is useful to comment on the differences between our approach to resummation 
and another method published recently by Catani et al (CMNT)§. We both use the 
same universal leading-logarithm expression in moment space, but differences occur 
after the transformation to momentum space. In this paper, we set aside comments on 
mathematical aspects of their procedure and focus instead on phenomenological issues 
of interest. As remarked above, the Mellin transformation generates sub leading terms 
in momentum space. CMNT choose to retain all of these inasmuch as they perform 
the Mellin inversion numerically. Instead, in keeping with the fact that subleading 
logarithmic terms are not universal, we retain only the universal leading logarithm 
terms in momentum space, and we restrict our phase space integration to the region 
in which the subleading terms would not be numerically significant regardless. The 
differences in the two approaches can be stated more explicitly if we examine the 
perturbative expansion of the kernel Hfj(z,a), Eq. ([[4]). If, instead of restricting 
the resummation to the universal leading logarithms only, we were to use the full 
content of Eq. (|T4l) , we would arrive at an analytic expression that is equivalent to 
the numerical inversion of CMNT, 



U% ~ 2aC ij 



£n 2 (l -z)- 2 lE tn(\ - z) 



0(a 2 ). (17) 



In terms of this expansion, in our work we retain only the leading term £n 2 (l — z) at 
order a, but CMNT retain both this term and the subleading term — 27#^n(l — z). 
We judge that it is not justified to keep the subleading term for three reasons: it is not 
universal; it is not the same as the subleading term in the exact 0(a 3 ) calculation; and 
it can be changed arbitrarily if one elects to keep non-leading terms in moment space. 
As a practical matter, the subleading term is negative and numerically significant. 
In the qq channel at m = 175 GeV and y/S = 1.8 TeV, its inclusion eliminates more 
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Fig. 2. Differential cross section da/dr) in the MS-scheme for the (a) and (b) gg channels: Born 
(dotted), NLO (dashed) and resummed (solid). 



than half of the contribution from the leading term. In our view, the presence of 
numerically significant subleading contributions means that non-universal structures 
are not under control and begs the question of consistency!!. A further justification 
for the retention of only the leading term is that it approximates the exact next-to- 
leading order result well, as shown in Fig. 1(a). The choice made by CMNT reproduces 
only one-third of the exact next-to- leading order result. The influence of subleading 
terms is amplified at higher orders where additional subleading structures occur in 
the CMNT approach with significant numerical coefficients proportional to 7r 2 , C(3), 
and so forth. Further comments about the different results in the two approaches are 
reserved to our discussion of predictions for cross sections at = 1.8 TeV. 



4. Physical cross section 



In order to achieve the best accuracy available we wish to include in our predictions 
as much as is known theoretically. Our final resummed partonic cross section can 
therefore be writteniH 



« pert 1 2 2 



a§ pert (r), m 2 , ft 2 ) - &§ +1 \r), m 2 , ft 2 ] 



+ a(° +1 W,/i 2 ) • (18) 



R;pert 



The second term is the part of the partonic cross section up to one-loop that is 
included in the resummation, while the last term is the exact one-loop cross sectionll. 
To obtain physical cross sections, we insert Eq. (|18D into Eq. ([]]), and we integrate over 
7]. Other than the top mass, the only undetermined scales are the QCD factorization 
and renormalization scales. We adopt a common value \x for both, and we vary this 
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Fig. 3. Inclusive cross section for top quark production in the MS-scheme. The dashed curves show 
our perturbative uncertainty band, while the solid curve is our central prediction: (a) = 1.8 TeV 
and (b) y/S = 2 TeV. 



scale over the interval \ijm G {0.5, 2} in order to evaluate the theoretical uncertainty 
of the numerical predictions. We use the CTEQ3M parton densitiesi. 



A quantity of phenomenological interest is the differential cross section 



d<Jij (S,rn 2 ,rf) 
dr) 



Its integral over r\ is the total cross section. In Fig. 2 we plot these distributions for 
m = 175 GeV, y/S = 1.8 TeV, and \i = m. The full range of rj extends to 25, 
but we display the behavior only in the near-threshold region where resummation is 
important. We observe that, at the energy of the Tevatron, resummation is significant 
for the qq channel and less so for the gg channel. In Fig. 1(a), the dotted curve shows 
that our final resummed cross section in the qq channel, after integration over all 
t], lies about half-way between the cross sections obtained from the near-threshold 
leading logarithms at orders 0(a 3 ) and 0(a A ). 

We show our total tt-pro duct ion cross section as a function of top mass in Fig. 3. 
The central value of our predictions is obtained with the choice /i/m = 1, and the 
lower and upper limits are the maximum and minimum of the cross section in the 
range of the hard scale fi/m 6 {0.5,2}. At m = 175 GeV, the full width of the 
uncertainty band is about 10% . We consider that the variation of the cross section 
over the range \ijm G {0.5,2} provides a good overall estimate of uncertainty. For 
comparison, we note that over the same range of fi, the strong coupling strength a 
varies by ±10% atm= 175 GeV. Our prediction of Fig. 3(a) is in agreement with the 
data0. We find a tl [m = 175 GeV, y/S = 1.8 TeV) = 5.52±g;^ pb. Our cross section 
is larger than the next-to- leading order value by about 9%. Using a different choice 
of parton densities^, we find a 4% difference in the central value of our prediction! 
at m — 175 GeV. A comparison of the predictions! in the MS and DIS factorization 
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schemes also shows a modest difference at the level of 4%. 

In Fig. 3(b) we present our predictions for an upgraded Tevatron operating at 
y/S = 2 TeV. We predict a tl {m = 175 GeV, VS = 2 TeV) = 7.56±g;Jg pb. The 2 pb 
increase in the predicted top quark cross section over its value at y/S =1.8 TeV is 
about a 37% gain. 

Two other groups have published predictions for the cross section. At m — 
175 GeV and y/s = 1.8 TeV, the three values are: ct**(BG§!) = 5.52^2 pb; 
a"(LSvN0) = 4.95±8;I8 pb; and <x"(CMNT@) = 4.75j#«j pb. From the purely nu- 
merical point of view, all three predictions agree within their estimates of theoreti- 
cal uncertainty. However, the resummation methods differ, as discussed above, the 
methods for estimating the uncertainties differ, and different parton sets are used. 
Comparing with LSvNa, we find that our central values are 10 — 14% larger, and 
our estimated theoretical uncertainty is 9 — 10% compared with their 28% — 20%. 
Our Born cross section, however, is about 3 — 5% larger than the LSvN Born cross 
section due to the different parton distributions used in the two calculations. Both 
the central value and the band of uncertainty of the LSvN predictions are sensitive to 
their undetermined infrared cutoffs. To estimate theoretical uncertainty, we use the 
standard yU-variation, whereas LSvN obtain their uncertainty primarily from varia- 
tions of their cutoffs. CMNT calculate a central value of the resummed cross section 
(also with fi/m = 1) that is less than 1% above the exact next-to-leading order value. 
As explained earlier, the suppression of the effects of resummation arises from the 
retention by CMNT of numerically significant non-universal subleading logarithmic 
terms. Indeed, if the subleading term — 27^£n(l — z) is discarded in Eq. (|T7D, the 
residuals 8ij/a^ LO defined by cmntI increase from 0.18% to 1.3% in the qq produc- 
tion channel and from 5.4% to 20.2% in the gg channelS After addition of the two 
channels, the total residual 6/a NLO grows from the negligible value of about 0.8% 
cited by CMNT to the value 3.5%. While still smaller than the increase of about 9% 
that we obtain, the increase of 3.5% vs. 0.8% shows the substantial influence of the 
subleading logarithmic terms retained by CMNT. We explain above the reasons that 
we judge that our method is preferable theoretically. 

Turning to pp scattering at the energies of the Large Hadron Collider (LHC) at 
CERN, we note a few significant differences from pp scattering at the energy of the 
Tevatron. The dominance of the qq production channel is replaced by gg dominance 
at the LHC. Owing to the much larger value of y/S, the near-threshold region in the 
subenergy variable is relatively less important, reducing the significance of initial- 
state soft gluon radiation. Lastly, physics in the region of large subenergy, where 
straightforward next-to-leading order QCD is also inadequate, becomes significant 
for ti production at LHC energies. Using the approach described in this paper, we 
estimate a tl {m = 175 GeV, VS = 14 TeV) = 760 pb. 
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5. Discussion and Conclusions 

In this paper, we summarize the calculation of the inclusive cross section for 
top quark production in perturbative QCD, including the resummation of initial- 
state gluon radiation to all orders. The advantages of the perturbative resummation 
method!!! are that there are no arbitrary infrared cutoffs, and there is a well-defined 
perturbative region of applicability where subleading logarithmic terms are numer- 
ically suppressed. Our theoretical analysis shows that perturbative resummation 
without a model for non-perturbative behavior is both possible and advantageous. In 
perturbative resummation, the perturbative region of phase space is separated cleanly 
from the region of non-perturbative behavior. The former is the region where large 
threshold corrections exponentiate but behave in a way that is perturbatively stable. 
The asymptotic character of the QCD perturbative series, including large multiplica- 
tive color factors, is flat, and excursions around the optimum number of perturbative 
terms does not create numerical instabilities or intolerable scale-dependence. Infrared 
renormalons are far away from the stability plateau and, even though their presence 
is essential for defining this plateau, they are of no numerical consequence in the 
perturbative regime. Large color factors, which are multiplicative in the exponent, 
enhance the infrared renormalon effects and contribute significantly to limiting the 
perturbative regime. 

Our resummed cross sections are about 9% above the next-to-leading order cross 
sections computed with the same parton distributions. The scale dependence of our 
cross section is fairly flat, resulting in a 9 — 10% theoretical uncertainty. This variation 
is smaller than the corresponding dependence of the next-to-leading cross section, as 
should be expected. In recent papers!, the authors state that the increase in cross 
section they find with their resummation method is of the order of 1% over next- 
to-leading order. The numerical difference in the two approaches boils down to the 
treatment of the subleading logarithms, which can easily shift the results by a few 
percent, if proper care is not taken. Our approach includes the universal leading log- 
arithms only while theirs includes non-universal subleading structures which produce 
the suppression they find. In our opinion, their treatment of the subleading structures 
is not correct. The issue is one to be decided by the theory community; it is not one 
for experimental resolution. 

Our theoretical analysis and the stability of our cross sections under [i variation 
provide confidence that our perturbative resummation procedure yields an accurate 
calculation of the inclusive top quark cross section at Tevatron energies and exhausts 
present understanding of the perturbative content of the theory. Our prediction agrees 
with data, within the large experimental uncertainties. 

Extending our calculation to much larger values of m than shown in Fig. 3, we 
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find that resummation in the principal qq channel produces enhancements over the 
next-to-leading order cross section of 21%, 26%, and 34%, respectively, at m — 500, 
600, and 700 GeV. The reason for the increase of the enhancements with mass at 
fixed energy is that the threshold region becomes increasingly dominant. Since the 
qq channel also dominates in the production of hadronic jets at very large values of 
transverse momenta, we suggest that on the order of 25% of the excess cross section 
reported by the CDF collaboration!^ may well be accounted for by resummation. 

This work was supported by the U.S. Department of Energy, Division of High 
Energy Physics, Contract No. W-31-109-ENG-38. 
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